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Abstract 

We introduce an algorithm design technique for a class of combinatorial optimization 
problems with concave costs. This technique yields a strongly polynomial primal-dual 
algorithm for a concave cost problem whenever such an algorithm exists for the fixed- 
charge counterpart of the problem. For many practical concave cost problems, the 
fixed-charge counterpart is a well-studied combinatorial optimization problem. Our 
technique preserves constant factor approximation ratios, as well as ratios that depend 
only on certain problem parameters, and exact algorithms yield exact algorithms. 

Using our technique, we obtain a new 1.61-approximation algorithm for the concave 
cost facility location problem. For inventory problems, we obtain a new exact algo- 
rithm for the economic lot-sizing problem with general concave ordering costs, and a 
4-approximation algorithm for the joint replenishment problem with general concave 
individual ordering costs. 

1 Introduction 

We introduce a general technique for designing strongly polynomial primal-dual algorithms 
for a class of combinatorial optimization problems with concave costs. We apply the tech- 
nique to study three such problems: the concave cost facility location problem, the economic 
lot-sizing problem with general concave ordering costs, and the joint replenishment problem 
with general concave individual ordering costs. 

In the second author's Ph.D. thesis [Str08] (see also [MS12]), we developed a general 
approach for approximating an optimization problem with a separable concave objective 
by an optimization problem with a piecewise-linear objective and the same feasible set. 
When we are minimizing a nonnegative cost function over a polyhedron, and would like 
the resulting problem to provide a 1 -|- e approximation to the original problem in optimal 
cost, the size of the resulting problem is polynomial in the size of the original problem 
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and linear in 1/e. This bound implies that a variety of polynomial-time exact algorithms, 
approximation algorithms, and polynomial-time heuristics for combinatorial optimization 
problems immediately yield fully polynomial-time approximation schemes, approximation 
algorithms, and polynomial-time heuristics for the corresponding concave cost problems. 

However, the piecewise-linear approach developed in [Str08] cannot fully address several 
difficulties involving concave cost combinatorial optimization problems. First, the approach 
adds a relative error of l-|-e in optimal cost. For example, using the approach together with 
an exact algorithm for the classical lot-sizing problem, we can obtain a fully polynomial- 
time approximation scheme for the lot-sizing problem with general concave ordering costs. 
However, there are exact algorithms for lot-sizing with general concave ordering costs [e.g. 
Wag60, AP93], making fully polynomial-time approximation schemes of limited interest. 

Second, suppose that we are computing near-optimal solutions to a concave cost problem 
by performing a 1 -|- e piecewise-linear approximation, and then using a heuristic for the 
resulting combinatorial optimization problem. We are facing a trade-off between choosing 
a larger value of e and introducing an additional approximation error, or choosing a smaller 
value of e and having to solve larger combinatorial optimization problems. For example, 
in [Str08] , we computed near-optimal solutions to large-scale concave cost multicommodity 
flow problems by performing piecewise-linear approximations with e = 1%, and then solving 
the resulting fixed-charge multicommodity fiow problems with a primal-dual heuristic. The 
primal-dual heuristic itself yielded an average approximation guarantee of 3.24%. Since we 
chose e = 1%, the overall approximation guarantee averaged 4.27%. If we were to choose 
e = 0.1% in an effort to lower the overall guarantee, the size of the resulting problems would 
increase by approximately a factor of 10. 

Third, in some cases, after we approximate the concave cost problem by a piecewise- 
linear problem, the resulting problem does not reduce polynomially to the corresponding 
combinatorial optimization problem. As a result, the piecewise-linear approach in [Str08] 
cannot obtain fully polynomial-time approximation schemes, approximation algorithms, 
and polynomial-time heuristics for the concave cost problem. For example, when we carry 
out a piecewise-linear approximation of the joint replenishment problem with general con- 
cave individual ordering costs, the resulting joint replenishment problem with piecewise- 
linear individual ordering costs can be reduced only to an exponentially-sized classical joint 
replenishment problem. 

These difficulties are inherent in any piecewise-linear approximation approach, and can- 
not be addressed fully without making use of the problem structure. 

The technique developed in this paper yields a strongly polynomial primal-dual algo- 
rithm for a concave cost problem whenever such an algorithm exists for the corresponding 
combinatorial optimization problem. The resulting algorithm runs directly on the concave 
cost problem, yet can be viewed as the original algorithm running on an exponentially or 
infinitely-sized combinatorial optimization problem. Therefore, exact algorithms yield ex- 
act algorithms, and constant factor approximation ratios arc preserved. Since the execution 
of the resulting algorithm mirrors that of the original algorithm, we can also expect the 
aposteriori approximation guarantees of heuristics to be similar in many cases. 
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1.1 Literature Review 



1.1.1 Concave Cost Facility Location 

In the classical facility location problem, there are m customers and n facihties. Each 
customer i has a demand di > 0, and needs to be connected to an open facihty to satisfy 
this demand. Connecting a customer z to a facihty j incurs a connection cost Cijdi; we 
assume that the connection costs are nonnegative and satisfy the metric inequality. Each 
facility j has an associated opening cost fj € M+. Let Xij = 1 if customer i is connected to 
facility j, and Xij = otherwise. Also let yj = 1 if facility j is open, and yj = otherwise. 
Then the total cost is Yl^=i fjVj + Z^i^i X^J=i CijdiXij- The goal is to assign each customer 
to one facility, while minimizing the total cost. 

The classical facility location problem is one of the fundamental problems in operations 
research [CNW90, NW99]. The reference book edited by Mirchandani and Francis [MF90] 
introduces and reviews the literature for a number of location problems, including classical 
facility location. Since in this paper, our main contributions to facility location problems 
are in the area of approximation algorithms, we next provide a brief survey of previous 
approximation algorithms for classical facility location. 

Hochbaum [Hoc82] showed that the greedy algorithm provides a O(logn) approxima- 
tion for this problem, even when the connection costs Cij are not metric. Shmoys et al. 
[STA97] gave the first constant-factor approximation algorithm, with a guarantee of 3.16. 
More recently, Jain et al. introduced primal-dual 1.861 and 1.61-approximation algorithms 
[JMM+03]. Sviridenko [Svi02] obtained a 1.582-approximation algorithm based on LP 
rounding. Mahdian et al [MYZ06] developed a 1.52-approximation algorithm that com- 
bines a primal-dual stage with a scaling stage. Currently, the best known ratio is 1.4991, 
achieved by an algorithm that employs a combination of LP rounding and primal-dual 
techniques, due to Byrka [Byr07]. 

Concerning complexity of approximation, the more general problem where the connec- 
tion costs need not be metric has the set cover problem as a special case, and therefore is not 
approximable to within a certain logarithmic factor unless P = NP [RS97]. The problem 
with metric costs does not have a polynomial-time approximation scheme unless P = NP, 
and is not approximable to within a factor of 1.463 unless NP C DTIME (r/'^(^°g'°s'7)) 
[GK99]. 

A central feature of location models is the economies of scale that can be achieved by 
connecting multiple customers to the same facility. The classical facility location problem 
models this effect by including a fixed charge fj for opening each facility j. As one of 
the simplest forms of concave functions, fixed charge costs enable the model to capture 
the trade-off between opening many facilities in order to decrease the connection costs 
and opening few facilities to decrease the facility costs. The concave cost facility location 
problem generalizes this model by assigning to each facility j a nondecreasing concave cost 
function : M_|_ — ^ R_|_, capturing a wider variety of phenomena than is possible with 
fixed charges. We assume without loss of generality that <^j(0) = for all j. The cost at 
facility j is a function of the total demand at j, that is diXij), and the total cost 

Researchers have studied the concave cost facility location problem since at least the 
1960's [KH63, FLR66]. Since it contains classical facility location as a special case, the 
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previously mentioned complexity results hold for this problem — the more general non-metric 
problem cannot be approximated to within a certain logarithmic factor unless P = NP, 
and the metric problem cannot be approximated to within a factor of 1.463 unless NP C 
DTIME (jjO{iogiogri)y 

To the best of our knowledge, previously the only constant factor approximation al- 
gorithm for concave cost facility location was obtained by Mahdian and Pal [MP03], who 
developed a 3 + e approximation algorithm based on local search. 

When the concave cost facility location problem has uniform demands, that is di = d2 = 
■ ■ ■ = djn, a wider variety of results become available. Hajiaghayi et al. [HMM03] obtained 
a 1.861-approximation algorithm. A number of results become available due to the fact 
that concave cost facility location with uniform demands can be reduced polynomially to 
classical facility location. For example, Hajiaghayi et al. [HMM03] and Mahdian et al. 
[MYZ06] described a 1.52-approximation algorithm. 

In the second author's Ph.D. thesis [Str08] (sec also [MS12]), wc obtain a 1.4991 -|- e 
approximation algorithm for concave cost facility location by using piecewise-linear approx- 
imation. The running time of this algorithm depends polynomially on 1/e; when e is fixed, 
the running time is not strongly polynomial. 

Independently, Romcijn et al. [RSSZIO] developed strongly polynomial 1.61 and 1.52- 
approximation algorithms for this problem, each with a running time of O(n'^logn). Here 
n is the higher of the number of customers and the number of facilities. They consider 
the algorithms for classical facility location from [JMM+03, MYZ06] through a greedy 
perspective. Since this paper uses a primal-diial perspective, establishing a connection 
between the research of Romeijn et al. and ours is an interesting question. 

1.1.2 Concave Cost Lot-Sizing 

In the classical lot-sizing problem, we have n discrete time periods, and a single item 
(sometimes referred to as a product, or commodity). In each time period t = l,...,n, 
there is a demand dt € M-(- for the product, and this demand must be supplied from 
product ordered at time t, or from product ordered at a time s < t and held until time t. 
In the inventory literature this requirement is known as no backlogging and no lost sales. 
The cost of placing an order at time t consists of a fixed cost ft £ M+ and a per-unit cost 
Ct G M+: ordering units costs ft + Q^j. Holding inventory from time t to time t + 1 
involves a per-unit holding cost ht € M+: holding ^t units costs ht$,t- The goal is to satisfy 
all demand, while minimizing the total ordering and holding cost. 

The classical lot-sizing problem is one of the basic problems in inventory management 
and was introduced by Manne [Man58], and Wagner and Whitin [WW58]. The literature 
on lot-sizing is extensive and here we provide only a brief survey of algorithmic results; for a 
broader overview, the reader may refer to the book by Pochet and Wolsey [P W06] . Wagner 
and Whitin [WW58] provided a O(ra^) algorithm under the assumption that q < ct-i+ht-i; 
this assumption is also known as the Wagner- Whitin condition, or the non-speculative 
condition. Zabcl [Zab64], and Eppen et al [EGP69] obtained O(n^) algorithms for the 
general case. Federgruen and Tzur [FT91], Wagelmans et al. [WvHK92], and Aggarwal 
and Park [AP93] independently obtained O(nlogn) algorithms for the general case. 

Krarup and Bilde [KB77] showed that integer programming formulation used in Section 
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3 is integral. Levi et al. [LRS06] also showed that this formulation is integral, and gave a 
primal-dual algorithm to compute an optimal solution. (They do not evaluate the running 
time of their algorithm.) 

The concave cost lot-sizing problem generalizes classical lot-sizing by replacing the fixed 
and per-unit ordering costs ft and Ct with nondecreasing concave cost functions 4>t : M_|_ 
M+. The cost of ordering units at time t is now (ptUt)- We assume without loss of 
generality that 0t(O) = for all t. This problem has also been studied since at least the 
1960's. Wagner [WagGO] obtained an exact algorithm for this problem. Aggarwal and Park 
[AP93] obtain another exact algorithm with a running time of O(n^). 

1.1.3 Concave Cost Joint Replenishment 

In the classical joint replenishment problem (JRP), we have n discrete time periods, and 
K items (which may also be referred to as products, or commodities). For each item fc, the 
set-up is similar to the classical lot-sizing problem. There is a demand G M+ of item k in 
time period t, and the demand must be satisfied from an order at time t, or from inventory 
held from orders at times before t. There is a per-unit cost € M_|- for holding a unit 
of item k from time t to t + 1. For each order of item k at time t, we incur a fixed cost 

G M_l_. Distinguishing the classical JRP from K separate classical lot-sizing problems 
is the fixed joint ordering cost — for each order at time we pay a fixed cost of /° G R+, 
independent of the number of items or units ordered at time t. Note that and do 
not depend on the time t. The goal is to satisfy all demand, while minimizing the total 
ordering and holding cost. 

The classical JRP is a basic model in inventory theory [Jon87, AE88]. The problem 
is NP-hard [AJR89]. When the number of items or number of time periods is fixed, the 
problem can be solved in polynomial time [Zan66, Vei69]. Federgruen and Tzur [FT94] 
developed a heuristic that computes 1 + e approximate solutions provided certain input 
parameters are bounded. Shen et al [SSLT] obtained a 0(logn -|- log i^T) approximation 
algorithm for the one-warehouse multi-retailer problem, which has the classical JRP as a 
special case. Levi ct al. [LRS06] provided the first constant factor approximation algo- 
rithm for the classical JRP, a 2-approximation primal-dual algorithm. Levi et al. [LRS05] 
obtained a 2.398-approximation algorithm for the one-warehouse multi-retailer problem. 
Levi and Sviridenko [LS06] improved the approximation guarantee for the one-warehouse 
multi-retailer problem to 1.8. 

The concave cost joint replenishment problem generalizes the classical JRP by re- 
placing the fixed individual ordering costs by nondecreasing concave cost functions 
0^ : R+ — > R+. We assume without loss of generality that 0*^(0) = for all k. The methods 
employed by Zangwill [Zan66] and Veinott [Vei69] for the classical JRP with a fixed number 
of items or fixed number of time periods can also be employed on the concave cost JRP. 
We are not aware of results for the concave cost JRP that go beyond those available for the 
classical JRP. Since prior to the work of Levi et al. [LRS06] , a constant factor approxima- 
tion algorithm for the classical JRP was not known, we conclude that no constant factor 
approximation algorithms are known for the concave cost JRP. 
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1.2 Our Contribution 



In Section 2, we develop our algorithm design technique for the concave cost facility location 
problem. In Section 2.1, we describe preliminary concepts. In Sections 2.2 and 2.3, we 
obtain the key technical insights on which our approach is based. In Section 2.4, we obtain 
a strongly polynomial 1.61-approximation algorithm for concave cost facility location with 
a running time of 0{m^n + mn log n). We can also obtain a strongly polynomial 1.861- 
approximation algorithm with a running time of 0{m^n + mn log n). Here m denotes the 
number of customers and n the number of facilities. 

In Section 3, we apply our technique to the concave cost lot-sizing problem. We first 
adapt the algorithm of Levi et al. [LRS06] to work for the classical lot-sizing problem as 
defined in this paper. Levi et al. derive their algorithm in a slightly different setting that 
is neither a generalization nor a special case of the setting in this paper. In Section 3.1, 
we obtain a strongly polynomial exact algorithm for concave cost lot-sizing with a running 
time of O(n^). Here n is the number of time periods. While the running time matches 
that of the fastest previous algorithm [AP93], our main goal is to use this algorithm as 
a stepping stone in the development of our approximation algorithm for the concave cost 
JRP in the following section. 

In Section 4, we apply our technique to the concave cost JRP. We first describe the 
difficulty in using piecewise-linear approximation on the concave cost JRP. We then intro- 
duce a more general version of the classical JRP, which we call generalized JRP, and an 
exponentially-sized integer programming formulation for it. In Section 4.1, using the 2- 
approximation algorithm of Levi et al. [LRS06] as the basis, we obtain an algorithm for the 
generalized JRP that provides a 4-approximation guarantee and has exponential running 
time. In Section 4.2, we obtain a strongly polynomial 4-approximation algorithm for the 
concave cost JRP. 



2 Concave Cost Facility Location 

We first develop our technique for concave cost facility location, and then apply it to other 
problems. We begin by describing the 1.61-approximation algorithm for classical facility 
location due to Jain et al. [JMM+03] . We assume the reader is familiar with the primal-dual 
method for approximation algorithms [see e.g. GW97]. 

Let [n] = {1, . . . The classical facility location problem, defined in Section 1.1.1, 
can be formulated as an integer program as follows: 

n m n 

j=l i=\ j=l 

n 

s.t. Y^Xij = 1, i E [m], (lb) 

< Xij < yj, ie[m\,j e [n], (Ic) 

%e{0,l}, jG[n]. (Id) 

Recall that fj G M_|_ are the facility opening costs, Cij G R-|_ are the costs of connecting 
customers to facilities, and di € M_|_ are the customer demands. We assume that the 
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connection costs Cij obey the metric inequality, and that the demands dj are positive. Note 
that we do not need the constrains Xij G {0, 1}, since for any fixed y G {0, 1}", the resulting 
feasible polyhedron is integral. 

Consider the linear programming relaxation of problem (1) obtained by replacing the 
constraints yj G {0, 1} with yj > 0. The dual of this LP relaxation is: 

(2a) 

i€[m],je[n], (2b) 

3 e H, (2c) 

ie[m\,j e [n]. (2d) 

Since Wij do not appear in the objective, we can assume that they are as small as possible 
without violating constraint (2b). In other words, we assume the invariant Wij = max{0, Vi — 
Cijdi}. We will refer to dual variable Vi as the budget of customer i. 11 Vi > Cijdi, we say 
that customer i contributes to facility j, and Wij is its contribution. The total contribution 
received by a facility j is YliLi "^ij- A. facility j is tight if YllLi "^ij — fj over-tight if 

Em ^ r 

The primal complementary slackness constraints are: 

Xij{vi - Cijdi - Wij) = 0, i e [m],j G [n], (3a) 

Suppose that {x,y) is an integral primal feasible solution, and {v,w) is a dual feasible 
solution. Then, constraint (3a) says that customer i can connect to facility j (i.e. Xij = 1) 
in the primal solution only if j is the closest to i with respect to the modified connection 
costs Cij + Wij/di- Constraint (3b) says that facility j can be opened in the primal solution 
(i.e. yj = 1) only if it is tight in the dual solution. 

The algorithm of Jain ct al. starts with dual feasible solution {v,w) = and itcrativcly 
updates it, while maintaining dual feasibility and increasing the dual objective. (The in- 
crease in the dual objective is not necessarily monotonic.) At the same time, guided by the 
primal complementary slackness constraints, the algorithm constructs an integral primal 
solution. The algorithm concludes when the integral primal solution becomes feasible; at 
this point the dual feasible solution provides a lower bound on the optimal value. 

We introduce the notion of time, and associate to each step of the algorithm the time 
when it occurs. In the algorithm, we denote the time by t. 



max 

i=l 



s.t. Vi < Cijdi + Wij, 



J2 '^io < fj 

i=l 

Wij > 0, 
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Algorithm FLPD(m, n€Z+;cG W^", f GWl,dG M!p) 

(1) Start at time i = with the dual sohition (f,u>) = 0. AU facilities are 
closed and all customers are unconnected, i.e. {x, y) = 0. 

(2) While there are unconnected customers: 

(3) Increase t continuously. At the same time increase Vi and Wij for un- 
connected customers i so as to maintain Vi = tdi and Wij = max{0, Vi — 
Cijdi}. The increase stops when a closed facility becomes tight, or an 
unconnected customer begins contributing to an open facility. 

(4) If a closed facility j became tight, open it. For each customer i that 
contributes to j, connect i to j, set vt = Cijdi, and set Wij' = max{0, Vi — 
Cij ' d{ } for all facilities j' . 

(5) If an unconnected customer i began contributing to an open facility j, 
connect i to j. 

(6) Return (a;, y) and 

In case of a tie between tight facilities in step (4), between customers in step (5), or 
between steps (4) and (5), we break the tie arbitrarily. Depending on the customers that 
remain unconnected, in the next iteration of loop (2), another one of the facilities involved 
in the tie may open immediately, or another one of the customers involved in the tie may 
connect immediately. 

Theorem 1 (JMM+03). Algorithm FLPD is a 1.61-approximation algorithm for the clas- 
sical facility location problem. 

Note that the integer program and the algorithm in our presentation are different from 
those in [JMM+03]. However both the integer program and the algorithm are equivalent 
to those in the original presentation. 

2.1 The Technique 

The concave cost facility location problem, also defined in Section 1.1.1, can be written as 
a mathematical program: 

n \ m n 

j +^^CijdiXij, (4a) 
/ i=i j=i 

i e[m], (4b) 

i € [m],j G [n]. (4c) 

Here, : R_|_ R_|. arc the facility cost functions, with each function being concave 
nondecreasing. Assume without loss of generality that 4>j{0) = for all cost functions. 
We omit the constraints Xij G {0, 1}, which are automatically satisfied at any vertex of 
the feasible polyhedron. Since the objective is concave, this problem always has a vertex 
optimal solution [HH61]. 

Suppose that the concave functions (pj are piecewise linear on (0, +oo) with P pieces. 





^ m 


j=i \ 




n 


= 1, 
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The functions can be written as 



\o, ^j = 0. ^ ' 

As is well-known [e.g. FLR66], in this case problem (4) can be written as the following 
integer program: 

n P m n P 

™™ fipVoP + ZlZmi ^^^i + Sjp)diXijp, (6a) 

j=i p=i 1=1 j=i p=i 

n P 

s.t. ^^Xjjp = 1, ie[m], (6b) 
j=i p=i 

< Xijp < yjp, ie[m\,j e[n],pe [P], (6c) 

yjpe {0,1}, je[n],pe[P]. (6d) 

This integer program is a classical facility location problem with Pn facilities and m cus- 
tomers. Every piece p in the cost function cpj of every facility j in problem (4) corresponds 
to a facility {j,p} in this problem. The new facility has opening cost fjp. The set of cus- 
tomers is the same, and the connection cost from facility {j,p} to customer i is Cij + Sjp. 
Note that the new connection costs again satisfy the metric inequality. 

We now return to the general case, when the functions (l>j need not be piccewise linear. 
Assume that 0^ are given by an oracle that returns the function value (pji^j) ^-iid derivative 
<p'jiS,j) in time 0(1) for > 0. If the derivative at does not exist, the oracle returns the 

right derivative, and we denote (pji^j) = lim^_^^^_)_ '^^ ^ ■ The right derivative always 
exists at > 0, since is concave on [0, -|-oo). 

We interpret each concave function as a picccwisc-linear function with an infinite 
number of pieces. For each p > 0, we introduce a tangent fjp + sjp^j to (pj at p, with 



Djp 



ip) , fjp = (p) - PSjp. (7) 



We also introduce a tangent fjo + SjoCj to (pj at 0, with fjo = limp_j.o+ fjp and sjo = 
limp_^o+ Sjp. The limit limp_>.o+ fjp is finite because fjp are nondecreasing in p and bounded 
from below. The limit limp_^o+ Sjp is either finite or -|-oo because Sjp are nonincreasing in 
p, and we assume that this limit is finite. 

Our technique also applies when limp_j.o+ Sjp = +oo, in which case we introduce tangents 
to (f)j only at points p > 0, and then proceed in similar fashion. In some computational 
settings, using derivatives is computationally expensive. In such cases, we can assume that 
the demands are rational, and let di = -4r with d'l > 0, and d'- and d'' coprime integers. 

i 

Also let A = ^„ . Then, we can use the quantity '^^^L^.'+^jj-'^-'^L^i-l) instead of ^^ (^j) 
throughout. 

The functions (j)j can now be expressed as: 

4>j{^j) = h^^f^^+'^^^^--p^^'^^ (8) 

' ' lO, ^j = 0. 
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When is linear on an interval [Ci, C2]5 all points p G [Ci, C2) yield the same tangent, that 
is {fjp,Sjp) = {fjq,Sjq) for any p,q £ [CiX2)- For convenience, we consider the tangents 
{fjp, Sjp) for all p > 0, regardless of the shape of (pj. Sometimes, we will refer to a tangent 
{fjp, Sjp) by the point p that gave rise to it. 

We apply formulation (6), and obtain a classical facility location problem with m cus- 
tomers and an infinite number of facilities. Each tangent p to cost function (pj of facility j 
in problem (4) corresponds to a facility {j, p} in the resulting problem. Due to their origin, 
we will sometimes refer to facilities in the resulting problem as tangents. 

The resulting integer program is: 



3=1 p>0 i=l 3=1 p>0 



S.t. Xijp — 1, 

j=l p>0 

^ Xijp ^ Vjpi 



i E [m\, 

i G [m\,j G [n\,p > 0, 
j G [n],p > 0. 



(9a) 

(9b) 

(9c) 
(9d) 



Of course, we cannot run Algorithm FLPD on this problem directly, as it is infinitcly-sizcd. 
Instead, we will show how to execute Algorithm FLPD on this problem implicitly. Formally, 
we will devise an algorithm that takes problem (4) as input, runs in polynomial time, and 
produces the same assignment of customers to facilities as if Algorithm FLPD were run on 
problem (9). Thereby, we will obtain a 1.61-approximation algorithm for problem (4). We 
will call the new algorithm ConcaveFLPD. 

The LP relaxation of problem (9) is obtained by replacing the constraints yjp G {0, 1} 
with yjp > 0. The dual of the LP relaxation is: 



max 



i=l 

s.t. Vi < {cij + Sjp)di + Wijp , 

m 



i=l 

Wijp > 0, 



i G [m], j G [n],p > 0, 
i G H,p>0, 
i G [m], j G [ra],p > 0. 



(10a) 
(10b) 
(10c) 
(lOd) 



Since the LP relaxation and its dual are infinitely-sized, the strong duality property does 
not hold automatically, as in the finite LP case. However, we do not need strong duality 
for our approach. We rely only on the fact that the optimal value of integer program (9) is 
at least that of its LP relaxation, and on weak duality between the LP relaxation and its 
dual. 



2.2 Analysis of a Single Facility 

In this section, we prove several key lemmas that will enable us to execute Algorithm 
FLPD implicitly. We prove the lemmas in a simplified setting when problem (4) has only 
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one facility, and all connection costs Cij are zero. To simplify the notation, we omit the 
facility subscript j. 

Imagine that we are at the beginning of step (3) of the algorithm. To execute this step, 
we need to compute the time when the increase in the dual variables stops. The increase 
may stop because a closed tangent became tight, or because an unconnected customer 
began contributing to an open tangent. We assume that there are no open tangents, which 
implies that the increase stops because a closed tangent became tight. 

Let t = at the beginning of the step, and imagine that t is increasing to +00. The 
customer budgets start at > and increase over time at rates Si. At time t, the budget 
for customer i has increased to Vi+tSi. Connected customers arc modeled by taking Si = 0, 
and unconnected customers by taking Si = di. Denote the set of connected customers by 
C and the set of unconnected customers by U, and let n = \U\. 

First, we consider the case when all customers have zero starting budgets. 

Lemma 1. If Vi = for i € [m], then tangent p* = X^jg^/dj becomes tight first, at time 

f • 

t* = Sp* + . // there is a tie, it is between at most two tangents. 
Proof. A given tangent p becomes tight at time Sp + , ■ Therefore, 



p = argmm 
p>0 



+ '^^ ^ j = argmin < Sp + /p > • (11) 



The quantity Sp X^jgjj di + fp can be viewed as the value of the affine function fp + Sp^ at 
^ = di- Since fp + Sp^ is tangent to (j), and (f) is concave, 

fp + Sp^di>(plj2'^i] fo"^ P^^- (12) 

On the other hand, for tangent p* = "^i^ijdi, we have fp* + Sp* J2ieu '^■i- ~ 'l^i^ieu ^i) ■ 
Therefore, tangent p* becomes tight first, at time t* = Sp* + (Sec Figure 1.) 

Concerning ties, for a tangent p to become tight first, it has to satisfy fp + Sp ^i^jj di = 
4>{J2i(zu di) , or in other words it has to be tangent to (f) at J2ieu consider two cases. 

First, let C2 be as large as possible so that (j) is linear on [p*, ^2]- Then, any point p G [p*, C2) 
yields the same tangent as p* , that is {fp,Sp) = {fp*,Sp*). Second, let (i be as small as 
possible so that 4> is linear on [Ci,P*]- Then, any point p € [Ci,P*) yields the same tangent 
as Ci; that is {fp,Sp) = (/^n-s^J. Tangent Ci is also tangent to 4> at YlieU^i^ 
be different from tangent p*. Tangents p ^ [Ci)C2] have fp + Sp^^^jj di > ^i^i^udi). 
Therefore, in a tie, at most two tangents, ^1 and p* , become tight first. □ 

Next, we return to the more general case when customers have nonnegative starting 
budgets. Define 

Pi{t) = min{p >Q:Vi + t5i> Spdi], i G [m], (13) 

If Vi + tSi < Spdi for every p > 0, let Pi{t) = +00. Otherwise, the minimum is well-defined, 
since Sp is right-continuous in p. 
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Figure 1: Illustration of the proof of Lemma 1. 



Intuitively, Pi{t) is the leftmost tangent to which customer i is contributing at time 
t. Note that Sp is decreasing in p, since i;^ is a concave function. Therefore, customer i 
contributes to every tangent to the right of Pi{t), and does not contribute to any tangent 
to the left pi{t). For any two customers i and j, 



{vi + t6i)/di > {vj + t6j)/dj =^ pi{t) < Pj{t), 
{vi + t6i)/di = {vj + t6j)/dj =^ pi{t) = Pj{t), 
{vi + t6i)/di < {vj + t6j)/dj pi{t) > pj{t). 



(14a) 
(14b) 
(14c) 



Assume without loss of generality that the set of customers is ordered so that customers 
1, . . . , fj, are unconnected, customers fi + 1, . . . ,m are connected, and 



vi/di > V2/d2 > • • • > v^/d^, 

V^+l/d^+l > Vf,+2/d^+2 > • • • > Vm/d„ 



(15a) 
(15b) 



Note that {vi + t5i)/di = vi/di for connected customers, and {vi + t6i)/di = Vi/di + t for 
unconnected ones. By property (14), at all times t, we have pi{t) < p2{t) < • • • < P^{t) and 
p^+i(t) < Pf^+2{t) < • • • < Pmi't)- As t increases, Pi{t) for i G C are unchanged, while Pi{t) 
for i G U decrease. (See Figure 2.) 
Let 



lUt) = [pkit),Pk+iit)), 

If{t) = [plit),pl+l{t)), 



1 <k < fi, 
H + 1 < I < m, 



(16a) 
(16b) 



with I^{t) = [0,pi(t)) and i;^it) = \pf,it),+oo), as well as r^{t) = [0,p^+i(t)) and /^(t) = 
\Pmit),+oo). When an interval has the form [+oo,+oo), we interpret it to be empty. 
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Figure 2: Illustration of the definition of Pi{t). Here U = {1,2} and C = {3,4}. The gray 
arrows show how Pi{t) change as t increases. The inequalities show the set of customers 
that contribute to the tangents in each of the intervals defined by Pi{t). 

Consider the intervals 

Iki{t) = lUt) n If{t), 0<k<fi<l<m. (17) 

At any given time t, some of the intervals Iki{t) may be empty. As time increases, these 
intervals may vary in size, empty intervals may become non-empty, and non-empty intervals 
may become empty. The intervals partition [0, -|-oo), that is ^Q<k<ij,<l<mlkl{t) = [0, +oo), 
and Iki{t) n Irs{t) = for {k,l) / (r,s). 

Let ujp{t) be the total contribution received by tangent p at time t. The tangents on each 
interval Iki{t) receive contributions from unconnected customers {1, . . . , /c} and connected 
customers + 1, . . . , /}. We define C{k,l) = {1, . . . ,k} Li {/jL + 1, . . . ,1} to be the set of 
customers that contribute to tangents in Iki{t). 

For each interval Iki{t) with A; > 1, we define an alternate setting A{k,l), where all 
starting budgets are zero, customers in C{k,l) increase their budgets at rates di, and the 
remaining customers do not change their budgets. Let OJp\Tki) be the total contribution 
received by tangent p at time Tki in the alternate setting A{k^ I). We establish a correspon- 
dence between times t in the original setting and times Tfc; in A{k, I), given by t^i = Pkl+ctklt: 
with aki = Ei=i di/ Y.i&c(k,i) di and /3m = EiGC(fc,0 ^i/ J2iGC{k,l) di- Since aki > 0, times 
t G [0, -|-oo) are mapped one-to-one to times t^i G [f^ki, +oo). 

The following two lemmas relate the original setting to the alternate settings A[k,l). 

Lemma 2. Given a time t and an interval I ki{t) with k >1, any tangent p G receives 
the same total contribution at time t in the original setting as at time tj-i in A{k^ I), that is 
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Proof. The total contribution to p at time t in the original setting is 

k I 

i=l i=ix+\ 
k 

ieC{k,l) i=l ieC{k,l) 

= {Pkl + oikit - Sp) ^ = {Tki - Sp) ^ di. 

ieC{k,l) ieC{k,l) 

Since cjp{t) > 0, it follows that t^i — Sp > 0, and therefore 

{Tkl-Sp) ^ di= ^ max{0, Tfeidj - SpcZj} =a;p'(Tjfci). (19) 
ieC{k,l) ieC{k,l) □ 

Lemma 3. Given a time t and an interval Iki{t) with k > 1, a tangent p receives at least 
as large a total contribution at time t in the original setting as at time t^i in A{k,l), that 
is Up{t) > uf{Tki). 

Proof. If Tu — Sp < 0, then cOp^Tki) = 0, and thus uip{t) > ujp^Tki). If Tki — Sp > 0, then 
^pim) = {Tki-Sp)Y,i^c{k,i)di = ULii^i + ^di- Spdi) + Yl\=^+iivi- Spdi). Letpe Irs{t) 
for some r and s, and note that ujp{t) = YH^ii'^^i + ~ Spdi) + Z^i=^+i(^^i — Spdi). 
The difference between the two contributions can be written as 

r k 

UJp{t) - Up {Tki) = X + ^'^^ ~ ^P'^') ~ X + ^'^^ ~ ^P'^') 
i=k+l i=r+l 

s I 

i=l+l i=s+l 

Note that at least two of the four summations in this expression are always empty. We now 
examine the summations one by one: 

1. X^i=fc+i(^t + ic?« — Spdi). This summation is nonempty when r > k. In this case, in the 
original setting, customers /c + 1, . . . ,r do contribute to tangents in Irs{t) at time t. 
Therefore, Vi + tdi — Spdi > for i = k + 1, . . . ,r, and the summation is nonnegative. 

2. — X]iLr+i(^« "I" ~ ^pdi)- This summation is nonempty when r < k. In this case, 
in the original setting, customers r + 1, . . . , A; do not contribute to tangents in Irs{t) 
at time t. Therefore, Vi + tdi — Spdi < for i = r + 1, . . . ,k, and the summation is 
nonnegative. 

3. '^i=i^i{vi — Spdi). This summation is nonempty when s > I. In this case, in the 
original setting, customers Z + 1, . . . , s do contribute to tangents in Irs{t) at time t. 
Therefore, Vi — Spdi > for i = Z + 1, . . . , s, and the summation is nonnegative. 

4. — 5^'=s_|_i(vi — Spdi). This summation is nonempty when s < I. In this case, in the 
original setting, customers s + 1, . . . ,1 do not contribute to tangents in /rs(i) at time 
t. Therefore, Vi — Spdi < for i = s + 1, . . . , Z, and the summation is nonnegative. 
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As a result of the above cases, we obtain that ojp{t) — Up^Tki) > 0. □ 

When A: > 1, we can apply Lemma 1 to compute the first tangent to become tight in 
A{k, I), and the time when this occurs. Denote the computed tangent and time by and 

t'^i, and note that p',^i = T.ieC{k,l) and r^^ = Sp^^ + Let = ^^^^ be the time in 
the original setting that corresponds to time rj^i in A{k,l). Let t* = minji^^ ■ ^ < k < n < 
I < m}, and p* = p' . ,., ^ i . The following two lemmas will enable us to show 
that tangent p* becomes tight first in the original setting, at time t* . 

Lemma 4. // o tangent p becomes tight at a time t in the original setting, then t >t* . 

Proof. Since UQ<r<fM<s<mirs{t) = [0, +oo), there is an interval Iki{t) that contains p. Since 
the contributions to tangents in the interval Ioi{t) do not increase over time, A; > 1. 

Tangent p is tight at time t in the original setting, and therefore C0p{t) = fp. By Lemma 

2, p £ Iki{t) implies that ujp^Tki) = fp, and hence p is tight at time Tki in A{k, I). It follows 
that Tki > T^;, and therefore t>t'f^i>t*. □ 

Lemma 5. Each tangent with k > 1 becomes tight at a time t < t'^^ in the original 
setting. 

Proof. Since tangent p'^i is tight at time Tj^^ in A{k, I), we have w*^/ (r^^) = /„/ . By Lemma 

Pkl 

3, ujpi^^{t'j^i) > fp'^^, which means that is tight or over-tight at time t'j^i in the original 
setting. Therefore, p'f.^ becomes tight at a time t < t'f^i in the original setting. □ 

We now obtain the main result of this section. 

Lemma 6. Tangent p* becomes tight first in the original setting, at time t* . The quantities 
p* and t* can be computed in time 0{m'^). 

Proof. Lemma 4 implies that tangents become tight only at times t > t* . Lemma 5 implies 
that tangent p* = Pargmin{t;^,,:i<fc</i<Km} becomes tight at a time t < min{i^^ ■ ^ < k < /j. < 
I < m} = t*. Therefore, tangent p* becomes tight first, at time t*. 

To evaluate the running time, note that the di and Vi can be sorted in 0{m log m) time. 
Once the di and Vi are sorted, we can compute all quantities aki and /3ki in 0{m?), and 
then compute all and p'j^^ in O(ra^) via Lemma 1. Therefore, the total running time is 
0(m2). □ 

In case of a tie, Lemma 6 enables us to compute one of the tangents that become tight 
first. It is possible to obtain additional results about ties, starting with that of Lemma 1. 
However, we do not need such results in this paper, as Algorithm FLPD, as well as the 
algorithms in Sections 3 and 4, allow us to break ties arbitrarily. 

For many primal-dual algorithms, we can perform the computation in Lemma 6 faster 
than in O(m^), by taking into account the details of how the algorithm increases the dual 
variables. We will illustrate this with three algorithms in Sections 2.4, 3, and 4. 
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2.3 Other Rules for Changing the Dual Variables 

In this section, we consider the same setting as in the previous one, but in addition allow 
each customer i to change its budget at an arbitrary rate Si > 0. The rate is no longer 
limited to the set {0,di}, and we assume that at least one customer has 5i > 0. The 
following results are not needed to obtain the algorithms in this paper. We include them 
since they embody a more general version of our approach, and may be useful in developing 
primal-dual algorithms in the future. 

Consider the quantities Pi{t) as defined in equation (13). Since 6i need not equal di, 
the order of the Pi{t) may change as t increases from to +oo. At any given time t, the 
Pi{t) divide [0, +oo) into at most m + 1 intervals. For each set of customers K C [m], we 
introduce an interval 

Ixit) = [aKit),bKit)) = max pi(t), mm pi(t)). (21) 

iieK i^K / 

If if = 0, we set aif(t) = 0, and '\i K = [m], we set 6i^(t) = +00. If ax(t) > hxify or 

o-Kify — bxit) = +00, we take Ixit) to be empty. Note that Uxc[m] ^^(0 = [0' +0°)) ^.nd 
^Kit) n lL{t) = $ ioi K ^ L. Any interval that is formed by the Pi{t) as t increases from 
to +00 is among the intervals Ixit)- The set of customers contributing to tangents on an 
interval Ixit) is precisely K. 

As in the previous section, for each interval Ixit) with "^i^x^i ^ define an al- 

ternate setting A{K), where all starting budgets are zero, customers in K increase their 
budgets at rates di, and the remaining customers keep their budgets unchanged. We de- 
note the total contribution received by tangent p at time tk in A{K) by ojp{TK)- The 
correspondence between times t in the original setting and times tk in A(K) is given by 

TK = f3K + (y.Kt, with UK = Y^ieK T^ieK 13k = Y^ieK ^^il Yji^K ^i- Since ax > 0, 

the correspondence is one-to-one between times t G [0, -|-oo) and tk G [/Skj+oo). 

Lemma 7. Given a time t and an interval Ixit) with Y2ieK^i ^' ionQC-nt p E Ik{^) 
receives the same total contribution at time t in the original setting as at time tk in A{K), 
that is ujp(t) = oj^{tk)- 

Proof. The total contribution in the original setting is 



= X]^^* ~ ^P^^^ = X]^'"' (^Ktdi - Spdi) 
ieK ieK 



(22) 



= {Pk + OiKt - Sp) '^di = {tk - Sp) '^di= tOp{TK). |-| 

ieK ieK 

Lemma 8. Given a time t and an interval lK{t) with ^i^^^i > 0, a tangent p receives at 
least as large a total contribution at time t in the original setting as at time tk in A{K), 
that is LOp{t) > ujp{TK)- 

Proof. If Tk — Sp < 0, then iOp{t) > uj^{tk)- If tk — Sp > 0, let p G lL{t) for some 
L C [m], and note that lOp{tk) = (tk - Sp)J2ieKdi = J2ieKi'"i + ~ ^v^i), while 



16 



The difference between the two contributions is 

^p{i) - (tk) = {vi + t6i - Spdi) - {vi + t6i - Spdi). (23) 

i£L\K i£K\L 

Since p G lL{t)^ in the original setting, customers in L contribute to tangent p at time t, 
and therefore Vi + t5i — Spdi > for z G L, which implies that jg/,\|^ (I'i + t6i — Spdi) > 0. 
Conversely, customers not in L do not contribute to p at time t, implying that Vi+tSi—Spdi < 
for z ^ L, and therefore J2ieK\Li''^i + t^i^pdi) < 0. As a result, u!p{t) — uip{TK) > 0. □ 

Unlike in the previous section, we have exponentially many alternative settings A{K). 
The following derivations will enable us to compute the first tangent to become tight in 
the original setting, and the time when this occurs using only a polynomial number of 
alternative settings. 

As t increases from to +oo, the order of the quantities {vi + t5i)/di may change. Since 
the quantities are linear in t, as t ^ +oo, they assume an order that no longer changes. We 
use this order to define a permutation 7r(+oo) = (7ri(+oo), . . . ,7rm(+oo)), with 7rj(+oo) = j 
meaning that {vj +t6j)dj is the i-th largest quantity. If two quantities are tied as t +oo, 
we break the tie arbitrarily. Similarly, for any time t G [0,+oo), we define a permutation 
Tr{t) = (7ri(t), . . . , TTmit))- In this case, if two quantities are tied, we break the tie according 
to 7r(+oo). For example, suppose that the two largest quantities at time t are tied, that 
they are {vi + t5i)/di = {v2 + i(^2)/c?2, and that 7ri(+oo) = 1 and 7rj(+c») = 2 with i < j. 
Then we take 7ri(t) = 1 and TT2{t) = 2. 

Compare two quantities 

{vi + t6i)/di vs. {vj + t5j)/dj. (24) 

If their order changes as t increases from to +oo, then there is a > such that the 
sign between the quantities is '<' on [0,^), '=' at 6, and '>' on (^, +oo), or vice-versa. 
Let 9i < ■ ■ ■ < Ojihe all such times when the sign between two quantities changes, and let 
^0 = and Or+i = +oo. Since there are m(m — l)/2 pairs of quantities, R < m{m — l)/2. 
The proof of the following lemma follows from these definitions. 

Lemma 9. As t increases from to +oo, the permutation 7r(t) changes at times 9i, . . . , Or. 
Moreover, 7T{t) is unchanged on the intervals [9r, Or+i) for r = 0, . . . , R. 

We now bound the number of intervals Ixit) that ever become nonempty. Let }C{t) = 
{{iTiit), Tri{t)} : i = 0, . . . , m} and /C = U^=o^(^r-)' and note that \IC{t)\ < m + 1 and 
|/C| < (m + l)(m(m - l)/2 + 1) = 0{m^). 

Lemma 10. As t increases from to +oo, only intervals Ixit) with K E K. ever become 
nonempty, that is {K : 3t >0 s.t. Ixit) 7^ 0} C A^. 

Proof. Fix a time t, and note that by property (14), we have P7ri(t)(^) ^ P-K2{t){'t) < • • • < 
PT^^{t){t). Therefore, the intervals Ixit) may be nonempty only when K G K,{t). Since 7r(t) 
is unchanged on the intervals [Or, 6'r+i) for r = 0, . . . , i?, if an interval Ixit) ever becomes 
nonempty, then K E )C. □ 
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As in the previous section, when J2ieK^i > 0' compute the first tangent to 

become tight in A{K), and the time when this occurs using Lemma 1. Let the computed 

tangent and time be = Yli^K '''k ~ ^p'^ + ^"^^ ^'k ~ '^'^aK^^ ^® time 

in the original setting corresponding to time rj^ in A{K). Next, we show that tangent p* = 
P'.r^rrAn{t'^:Y.,eK5i>^,KeK} "^^^^^^^ ^^^"^^ ^^st, at time t* = min{i'^ : Y^i^K^i > 0,K G /C}. 

Lemma 11. If a tangent p becomes tight at a time t in the original setting, then t>t*. 

Proof. Let Ixit) be the interval that contains p. Since this interval is nonempty, K G IC, 
and since the contribution to p must be increasing over time, J2ieK > ^■ 

Tangent p is tight at time t in the original setting, and therefore ojp{t) = fp. By Lemma 
7, u!p (tk) = fp, and hence p is tight at time tk in A{K). It follows that tk > t^^, and 
therefore t>t'j^>t*. □ 

Lemma 12. Each tangent p'j^ with Y^^^ > ^ becomes tight at a time t < t'j^ in the 
original setting. 

Proof. Since p'j^ is tight at time in A{K), we have cu^ (rj^) = fp'^. By Lemma 8, 

C0p'^{t'^) > fp'^, which means that p'^ is tight or over-tight at time in the original 
setting. Therefore, becomes tight at a time t <t'j^ in the original setting. □ 

Lemma 13. Tangent p* becomes tight first in the original setting, at time t*. The quantities 
p* and t* can be computed in time O(m^). 

Proof. By Lemma 11, tangents only become tight at times t > t*, while by Lemma 12, p* = 
p' . r,, V- r ^ v-T becomes tight at a time t < minjt'^ : JZicr^Si > 0,K E K,} = t*. 

Therefore, p* becomes tight first, at time t*. 

Concerning the running time, note that the Or can be computed in 0{m?) and sorted in 
0{m? log m) time. The permutations 7r(+oo) and 7r(0) can be computed in 0{m log m) time. 
Processing the Or in increasing order, we can compute each 7r(^r) in 0{m) time amortized 
over all Or- Computing and for all K e lC{Or) takes 0{m) time. Therefore, the total 
running time is O(m^). □ 

The results in this section can be generalized further to allow the rates Si to be negative. 



2.4 Analysis of Multiple Facilities 

Wc now show how to execute Algorithm FLPD implicitly when problem (4) has multiple 
facilities. In Section 2.2, in addition to assuming the presence of only one facility, we 
assumed that the connection costs Cij were 0. We remove this assumptions as well. 

In this section, we continue to refer to facilities of infinitely-sized problem (9) as tan- 
gents, and reserve the term facility for facilities of concave cost problem (4). We say that 
customer i contributes to concave cost facility j if Vi > Cij. Wc distinguish between when 
a customer contributes to a concave cost facility j and when the customer contributes to a 
tangent p belonging to concave cost facility j. 
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When executing Algorithm FLPD imphcitly, the input consists of m, n, the connection 
costs Cjj, the demands d-i, and the cost functions (f)j, given by an oracle. As intermediate 
variables, we maintain the time t, and the vectors v, x, and y. For x and y, we maintain 
only the non-zero entries. The algorithm returns v, x, and y. We also maintain standard 
data structures to manipulate these quantities as necessary. Note that we do not maintain 
nor return the vector w, as any one of its entries can be computed through the invariant 
Wijp = max{0,Wj - {cij + Sjp)di}. 

Clearly, step (1) can be executed in polynomial time. In order to use induction, suppose 
that we have executed at most m — 1 iterations of loop (2) so far. Since the algorithm opens 
at most one tangent at each iteration, at any point at most m — 1 tangents are open. To 
analyze step (3) , wc consider three events that may occur as this step is executed: 

1. A closed tangent becomes tight. 

2. An unconnected customer begins contributing to an open tangent. 

3. An unconnected customer begins contributing to a facility. 

When step (3) is executed, the time t stops increasing when event 1 or 2 occurs. For 
the purpose of analyzing this step, we assume that t increases to +00 and that Vi for 
unconnected customers are increased so as to maintain Vi = tdi. 

Lemma 14. Suppose that event e at facility j is the first to occur after the beginning of 
step (3). Then we can compute the time t' when this event occurs in polynomial time. 

Proof. If e = 1, we use Lemma 6 to compute t'. The lemma's assumptions can be satisfied 
as follows. Since no events occur at other facilities until time t' , wc can assume that j is 
the only facility. Since the set of customers contributing to facility j will not change until 
time t', we can satisfy the assumption that Cij = by subtracting Cij from each Vi having 
Vi > Cij. Since an unconnected customer will not begin contributing to an open tangent 
until time t', we can assume that there are no open tangents. We can satisfy the assumption 
that t = at the beginning of step (3) by adding tdi to each Vi. 

If e = 2, wc compute t' by iterating over all unconnected customers and open tangents 
of facility j. If e = 3, we compute t' by iterating over all unconnected customers. □ 

When other events occur between the beginning of step (3) and time t', the computation 
in this lemma may be incorrect, however we can still perform it. Let t'^{j) be the time 
computed in this manner for a given e and j, and let t* = min{tg(j) : e G [3],j G [n]} and 
{e*,j*) = argmin{i;(j) : e G [S\,j G [n]}. 

Lemma 15. Event e* at facility j* is the first to occur after the beginning of step (3) . This 
event occurs at time t* . 

Proof. Suppose that an event e' at a facility j' occurs at a time t' <t*. If e' G {2, 3}, then 
t' > t'^,[j') > t* . This is a contradiction, and therefore this case cannot occur. 

If e' = 1, then we consider two cases. If there is an event e" G {2,3} that occurs at a 
facility j" at a time t" < t', then we use t" > t'^„{j") > t* to obtain a contradiction. If 
there is no such event e", then no new customer begins contributing to facility j' between 
the beginning of step (3) and time t'. Therefore, t' > t'^,{j') > t*, and we again obtain a 
contradiction. □ 
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Once we have computed t*, e*, and j*, we finish executing step (3) as fohows. If e* = 3, 
that is if the first event to occur is an unconnected customer beginning to contribute to 
J*, we update the hst of customers contributing to j* and recompute t* , e* , and j* . Since 
there are n facilities and at most m unconnected customers, event 3 can occur at most mn 
times before event 1 or 2 takes place. 

Once event 1 or 2 takes place, step (3) is complete, and we have to execute step (4) 
or (5). It is easy to see that these steps can be executed in polynomial time. Therefore, 
an additional iteration of loop (2) can be executed in polynomial time. By induction, 
each of the first m iterations of loop (2) can be executed in polynomial time. At each 
iteration, an unconnected customer is connected, either in step (4) or (5). Therefore, loop 
(2) iterates at most m times. Obviously, step (6) can be executed in polynomial time, and 
therefore Algorithm FLPD can be executed implicitly in polynomial time. Recall that we 
called the algorithm obtained by executing FLPD implicitly on infinitely-sized problem (9) 
ConcaveFLPD. 

Theorem 2. Algorithm ConcaveFLPD is a 1.61- approximation algorithm for concave 
cost facility location, with a running time of 0{m^n + mn log n). 

Proof. At the beginning of the algorithm, we sort the connection costs Cij, which can be 
done in 0{mn\og{mn)) time. Next, we bound the time needed for one iteration of loop 
(2). Note that since loop (2) iterates at most m times, there are at most m open tangents 
at any point in the algorithm. 

In step (3), we first compute min{t'j^(j) : j £ [n]}. Computing each t[{j) requires 0(rn?) 
per facility, and thus this part takes 0{m^n) overall. Next, we compute min{t2(i) : j £ [«]}, 
using 0(1) per customer and open tangent, and thus 0{vn?') overall. Finally, we compute 
vai\i{t'^{j) : j G [n]}. Since m.\ii{t'^[j) : j S [n]} = minjcij : Cij > t}, we have sorted the 
values Cij, and t only increases as the algorithm runs, this operation takes 0{mn) over 
the entire run of the algorithm. Therefore, we can determine the next event to occur in 
O(m^n). 

If event 1 or 2 is the next one, step (3) is complete. If event 3 is next, an uncon- 
nected customer begins to contribute to facility j* . In this case, we recompute and 
min{t'^(j) : j G [n]}. Recomputing can be done in 0(m), since we only have to add 

one customer to the setting of Lemma 6. Recomputing m.m.{t'i{j) : j G [n]} takes 0(1), as 
t'lij*) does not increase when an unconnected customer begins contributing to facility j* . 
Note that min{t2(j) : j G [n]} does not change. Next, we recompute min{t3(j) : j G [n]}, 
which takes 0{mn) over the entire run of the algorithm. The total time to process event 
3 and determine the next event to occur is 0(m). Event 3 occurs at most mn times be- 
fore event 1 or 2 occurs, and therefore the total time for processing event 3 occurrences is 
O(m^ra). 

Step (4) can be done in 0(m), and step (5) in 0(1). Therefore, the time for one iteration 
of loop (2) is O(m^n). Since there are at most m iterations of loop (2), the running time 
of the algorithm is 0{m^n + mn log n). 

By Theorem 1, Algorithm FLPD is a 1.61-approximation algorithm for problem (1). 
The approximation ratio for problem (4) follows directly from the fact that we execute 
Algorithm FLPD implicitly on infinitely-sized problem (9). □ 
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By a similar application of our technique to the 1.861-approximation algorithm for 
classical facility location of Mahdian et al. [JMM"'"03], we obtain a 1.861-approximation 
algorithm for concave cost facility location with a running time of 0{m^n + mnlogn). 

3 Concave Cost Lot-Sizing 

In this section, we apply the technique developed in Section 2 to concave cost lot-sizing. 
The classical lot-sizing problem is defined in Section 1.1.2, and can be written as a linear 
program: 

n n n 

min ^ fsVs + ^'^ics + hst)dtXst, (25a) 

s=l s=l t=s 

t 

s.t.^Xst = l, l<t<n, (25b) 

s=l 

0<Xst<ys, l<s<t<n. (25c) 

Recall that ft G M_|_ and Ct G M+ are the fixed and per-unit costs of placing an order at time 
t, and dt G M+ is the demand at time t. The per-unit holding cost at time t is /it G M+, and 
for convenience, we defined hgf = ^*Z^ hi. Note that we omit the constraints G {0, 1}, 
as there is always an optimal extreme point solution that satisfies them [KB77]. 

We now adapt the algorithm of Levi et al. [LRS06] to work in the setting of problem 
(25). Levi et al. derive their algorithm in a slightly different setting, where the costs hgt 
are not necessarily the sum of period holding costs ht, but rather satisfy an additional 
monotonicity condition. 

The dual of problem (25) is given by: 

(26a) 

l<s<t<n, (26b) 
1 < s < n, (26c) 

l<s<t<n. (26d) 

As with facility location, since the variables Wgt do not appear in the objective, we assume 
the invariant Wgt = max{0, vt — {cg + hst)dt}. Note that lot-sizing orders correspond to facil- 
ities in the facility location problem, and lot-sizing demand points correspond to customers 
in the facility location problem. 

We refer to dual variable vt as the budget of demand point t. If Vt > (cg + hst)dt, 
we say that demand point t contributes to order s, and Wgt is its contribution. The total 
contribution received by an order s is ^t=s^st- An order t is tight if Y17=s''^st = fs and 
over-tight if X^"^^ Wgt > fs- 



n 



max Vf, 

t=i 

s.t. Vt < {cs + hst)dt + Wst, 

n 

XI "^st < fs, 



t=s 

Wst > 0, 
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The primal complementary slackness constraints are: 

Xst{vt - (cs + hst)dt - Wst) = 0, I < s <t <n, (27a) 

Vs Wst -fs^=0, l<s<n. (27b) 

Let (x, y) be an integral primal feasible solution, and (v, w) be a dual feasible solution. 
Constraint (27a) says that demand point t can be served from order s in the primal solution 
only if s is the closest to t with respect to the modified costs Cg + hst + Wst/dt- Constraint 
(27b) says that order t can be placed in the primal solution only if it is tight in the dual 
solution. 

The algorithm of Levi et al., as adapted here, starts with dual feasible solution (w, w) = 
and iteratively updates it, while maintaining dual feasibility and increasing the dual 
objective. At the same time, guided by the primal complementary slackness constraints, 
the algorithm constructs an integral primal solution. The algorithm concludes when the 
integral primal solution becomes feasible. An additional postprocessing step decreases the 
cost of the primal solution to the point where it equals that of the dual solution. At this 
point, the algorithm has computed an optimal solution to the lot-sizing problem. 

We introduce the notion of a wave, which corresponds to the notion of time in the primal- 
dual algorithm for facility location. In the algorithm, we will denote the wave position by 
W , and it will decrease continuously from hin to 0, and then possibly to a negative value 
not less than — ci — f\. We associate to each step of the algorithm the wave position when 
it occurred. 

Algorithm LSPD(n & 'L+] c, f,d & M.X,h e K+"^) 

(1) Start with the wave sA, W = hin and the dual solution {v,w) = 0. All 
orders are closed, and all demand points are unserved, i.e. {x, y) = 0. 

(2) While there are unserved demand points: 

(3) Decrease W continuously. At the same time increase vt and Wst for 
unserved demand points t so as to maintain vt = max{0, dtihu — W)} 
and Wst = max{0, U4 — (cg -I- hat)dt}. The wave stops when an order 
becomes tight. 

(4) Open the order s that became tight. For each unserved demand point 
t contributing to s, serve t from s. 

(5) For each open order .s from 1 to n; 

(6) If there is a demand point t that contributes to s and to another open 
order s' with s' < s, close s. Reassign all demand points previously 

served from s to s' . 

(7) Return {x,y) and {v,w). 

In case of a tie between order points in step (4) , we break the tie arbitrarily. Depending 
on the demand points that remain unserved, another one of the tied orders may open 
immediately in the next iteration of loop (2). 

The proof of the following theorem is almost identical to that from [LRS06] , and there- 
fore for this proof we assume the reader is familiar with the lot-sizing results from [LRS06] . 

Theorem 3. Algorithm LSPD is an exact algorithm for the classical lot-sizing problem. 
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Proof. We will show that after we have considered open order s, at the end of step (6), we 
maintain two invariants. First, each demand point is contributing to the fixed cost of at 
most one open order from the set {1, . . . , s}. Second, each demand point is assigned to an 
open order and contributes to its fixed cost. 

The first invariant follows from the definition of the algorithm. Indeed, if a demand 
point t' is contributing to s' and s with s' < s, then the algorithm would have closed s. 

Clearly the second invariant holds at the beginning of loop (5). It continues to hold 
after we review order s if we have not closed s. Let us now consider the case when we 
have closed s. The demand points that have contributed to s can be classified into two 
categories. The first category contains the demand points whose dual variables stopped 
due to s becoming tight — these demand points were served from s and are now served from 
s'. Since t contributes to s', so do these demand points. The second category contains the 
demand points whose dual variables stopped due to another order s" becoming tight. The 
case s" < s cannot happen, or s would have never opened. Hence, s < s", and therefore 
s" is currently open. Moreover, these demand points are currently served from s" and are 
contributing to it. 

Therefore, at the end of loop (5), each demand point is contributing to the fixed cost of 
at most one open order. Therefore, the fixed cost of opening orders is fully paid for by the 
dual solution. Moreover, each demand point is served from an open order, and therefore 
the primal solution is feasible. Since each demand point contributes to the fixed cost of the 
order it is served from, the holding and variable connection cost is also fully paid for by 
the dual solution. Since the primal and dual solutions have the same cost, the algorithm is 
exact. □ 

3.1 Applying the Technique 

We now proceed to develop an exact primal-dual algorithm for concave cost lot-sizing. The 
concave cost lot-sizing problem is defined in Section 1.1.2: 

n / n \ n n 

min XI '^M XI J + X X ^stdtXst, (28a) 

s=l \t=s / s=l t=s 

t 

S.' 

s=l 

Xst>0, l<s<t<n. (28c) 



i.t.'^Xst = l, l<t<n, (28b) 



Here, the cost of placing an order at time t is given by a nondecreasing concave cost function 
(f>t : — > R+. We assume without loss of generality that (j)t{0) = for all t. 

The application of our technique to the lot-sizing problem is similar to its application to 
the facility location problem in Section 2. First, we reduce concave cost lot-sizing problem 
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(28) to the following infinitely-sized classical lot-sizing problem. 



n n n 



mm 



(29a) 



s=l p>0 s=l t=s p>0 



S.t. ^ ^ ^ ^ Xgpf — 



1, 



1 < t < n, 



(29b) 



^ ^^spt ^ ysp) 



1 < s < t < n,p > 0. 



(29c) 



Again we note that since LP (29) is infinitely-sized, strongly duality does not hold auto- 
matically for it and its dual. However, the proof of Algorithm LSPD relies only on weak 
duality. The fact that the algorithm produces a primal solution and a dual solution with 
the same cost implies that both solutions are optimal and that strong duality holds. 

Following Section 2, let ConcaveLSPD be the algorithm obtained by executing Algo- 
rithm LSPD implicitly on infinitely-sized problem (29). 

Theorem 4. Algorithm ConcaveLSPD is an exact algorithm for concave cost lot-sizing, 
with a running time ofO{n'^). 

Proof. We consider the following events that may occur as step (3) of Algorithm LSPD is 
executed: 



If for an order point t, no tangents become tight in the course of the algorithm, wc let 
W2{t) = +00. The wave positions Wi{t) can be computed for all t at the beginning of the 
algorithm in 0{n). 

We compute the positions W2{t) by employing a set of intermediate values W2{t). Each 
value 14^2(0 defined as the time when a tangent of order point t becomes tight in a 
truncated problem consisting of time periods i, t + 1, . . . , n. We compute a subset of these 
values as follows. First, we compute W2{n), which requires 0(1) time by Lemma 1. To 
compute W2{t) given that W2{t -I- 1), ... , (^) computed, we can employ Lemma 6. 

The dual variables representing demand points t, . . . ,n can be divided into three con- 
secutive intervals. First are the dual variables that are increasing at the same rate as part 
of the wave, then the dual variables Vk that are not increasing but exceed htk, and finally 
the dual variables vj- that are not increasing, do not exceed htk, and therefore play no role 
in this computation. We employ Lemma 6 and distinguish two cases: 

1. Lemma 6 can be used to detect if a tangent is overtight. This indicates that ^^^(O 
an earlier wave position than W2{t -|- 1), . . . , W^l^) some k. In this case, we delete 
W2{t + 1) from our subset and repeat the computation of W^IO as if order point t+1 
does not exist. 

2. There are no overtight tangents. Thus, a tangent becomes tight at a wave position 
less than or equal to 14^2 + !)• 1^ this case we set W2{t) to this wave position, and 
proceed to the computation of 1^2 ~ 



Time 



Ea<-ii1 



W2{t) 



The wave reaches demand point t, i.e. W — hit- 
A tangent p of order point t becomes tight. 
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After computing W2{t), consider the values that remain in our subset and denote them 
by W2{t), W^(7r(l)), . . . , VF^(7r(A;)) for some k. By induction, these values yield the correct 
times when tangents become tight for the truncated problem consisting of time periods 
t, . . . ,n. After we have computed VF2(1), the values VF2(t) remaining in our subset yield 
the correct times W2{t), with the other values W2{t) = +00. Therefore, loop (2) is complete. 

A computation by Lemma 6 requires O(n^) time in the worst case. Since in this setting, 
all dual variables that are increasing exceed all dual variables that are stopped, each ^^(t) 
can be computed by Lemma 6 in 0(n). Each time we use Lemma 6 for a computation, 
a value W2{t) is either removed from the list or inserted into the list. Since each value is 
inserted into the list only once, the total number of computations is 0(n), and the total 
running time for loop (2) is 0{v?). 

At the beginning of step (5), there are at most n open tangents, and n demand points, 
and therefore this loop can be implemented in O(n^) as well. □ 

Note that the values W2{t) also yield a dual optimal solution to the infinitely-sized LP. 
The solution can be computed from the W2{t)-s in time 0(n) by taking vt = hu — W2{(T{t)), 
where a{t) is the latest time period less than or equal to t that has W2{(T{t)) < +00. 



4 Concave Cost Joint Replenishment 

In this section, we apply our technique to the concave cost joint replenishment problem 
(JRP). The classical JRP is defined in Section 1.1.3, and can be formulated as an integer 
program: 

n n K n n K 

-in Y: fVs + E E f'ys + E E E ^U^^'t, (30a) 

s=l s=l k=l s=l t=s k=l 

t 

s.t. E = 1' l<t<n,k£ [K], (30b) 

s=l 

< a;Jt < 2/°, 1 < s <t <n,k e[K], (30c) 

Q<x%<y'^, l<s<t<n,ke[K], (30d) 

t/° G {0, 1}, y,^ G {0, 1}, l<s<n,ke [K] . (30e) 

Recall that /° € M+ is the fixed joint ordering cost, /^^ G M+ is the fixed individual ordering 
cost for item k, and G M+ is the demand for item k at time t. The per-unit holding cost 
for item k at time t is h^; for convenience we defined hg^ = X]i=s ^^i- 

The concave cost JRP, also defined in Section 1.1.3, can be written as a mathematical 
program as follows: 

n / n K \ n K / n \ n n K 

E E E d^-st + E E E + E E E ^'t^-st, m^) 



s=l \t=s k=l / s=l k=l \t=s / s=l t=s k=l 

t 

s ■ 

s=l 



t. ^x^^ = l, l<t<n,k&[K], (31b) 

s=l 

x^i > 0, l<s<t<n,ke [K]. (31c) 



25 



Here the individual ordering cost for item k at time t is given by a nondecreasing concave 
function cj)^ : Wc assume without loss of generality that (/)''{0) = for all k. 

The joint ordering cost at time t is given by the function (p^ : IR+ — t- M+. To reflect the fact 
that only the individual ordering costs are general concave, (p^ has the form ^''(0) = and 
00(0 = /° for e > 0. 

Consider the case when the individual ordering cost functions (f)^ are piecewise linear 
with P pieces: 



0) — 0' 



Unlike with concave cost facility location and concave cost lot-sizing, the piecewise-linear 
concave cost JRP does not reduce polynomially to the classical JRP. Since there are multiple 
items, different pieces of the individual ordering cost functions may be employed by 
different items k as part of the same order at time t. When each cost function consists of 
P pieces, we would need P^ time periods to represent each possible combination, thereby 
leading to an exponentially-sized IP formulation. 

We could devise a polynomially-sized IP formulation for the piecewise-linear concave 
cost JRP, however such a formulation would have a different structure from the classical 
JRP, and would not enable us to apply our technique together with the primal-dual algo- 
rithm of Levi et al. [LRS06] for the classical JRP. Instead, we reduce the piecewise-linear 
concave cost JRP to the following exponentially-sized integer programming formulation, 
which we call the generalized joint replenishment problem. Let tt = (pi, . . . ,Pk), and let 
[P]"" = {{pi,...,PK):Pie[P]}. 



mm 



E E fLys.+ E + (33a) 

se[n] se[n],k£[K] l<s<t<n 

7rG[Pp 7rG[P]^ fee[is:],7re[P]^ 

s.t. E "^snt = 1> l<t<n,ke[K], (33b) 

se[t] 

< < yl, l<s<t<n,kG[K],7T G [Pf, (33c) 

< xtt <yL l<s<t<n,ke[K],TT G [Pf, (33d) 

G {0, 1}, 2/,^^ G {0, 1}, l<s<n,ke[K],7re [Pf. (33e) 

The intuition behind the generalized JRP is that each time period t in the piecewise-linear 
concave cost JRP corresponds to P^ time periods {t, vr) in the generalized JRP. Each time 
period {t,ir) allows us to use a different combination tt = (pi, ■ ■ ■ ,Pk) of pieces of the 
individual order cost functions (p^, . . . , (f)^ . 

This formulation docs not satisfy the cost assumptions required for the 2-approximation 
algorithm of Levi et al. [LRS06]. In the next section, we will devise, starting from the algo- 
rithm of Levi et al, an algorithm for the generalized JRP that provides a 4-approximation 
guarantee and runs in exponential time. In Section 4.2, we will employ our technique to 
obtain a strongly polynomial 4-approximation algorithm for the concave cost JRP. 
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4.1 An Algorithm for the Generahzed JRP 

Consider the LP relaxation of IP (33) obtained by replacing the constraints G {0, 1}, 
G {0, 1} with > O,?/^^^ ^ 0- The dual of this LP relaxation is: 

K n 

max (34a) 

k=l t=l 

s.t. < (4 + h'gM^ + + nL, '-ll^i^!''''' (34b) 



< f'^' l<s<n,fce[ii-], . . X 

t=s 

K n 

< /°, 1 < . < n, TT G [P]^, (34d) 

fe=l t=s 

> n 7/^= > n l<s<t<n,k€[K], . X 



Since now both w^^^ and Ug^^ are not present in the objective, the invariants for them 



become more involved. When Yl^=s max{0, — {cll + h'^t)(^t} ^ /pi,) '^^ let as before 



w'^^t = max{0, v"^ - (4 + 4)4} < 4^^. (35a) 

When X]r=s max{0, — (c^^ + h%)d^} > fp^, the algorithm will have fixed the values Wg^^ 
at the point when J2t=s max{0, — (c^^ + hgf)df} = f^^. In this situation, we let 

«L = max{0, i;*^ - (^ + /i.^J^ - «;,^^J. (35b) 

We now have demand points for every time-item pair, and we refer to as the budget of 
item k at time t. Given tt, if > {cp^ + hg^)dt, we say that demand point [k, t) contributes 
to the fixed cost of individual order {s,k,n) and Wg^f is its contribution. If > (c^^ + hfltjdt 
and X^tLs ^sTTt = /pfe' '^^ t^^t demand point (A;, t) contributes to the fixed cost of joint 
order (s, tt) and nj^^ is its contribution. 

Since we now have several items, each with its own holding costs, we think of as a 
"master" wave, and decrease it from n to 1 and then to a bounded amount below 1. For 
each item k, we maintain an item wave 

= h^\w\+h\_w\{W -[W\). (36) 

Intuitively, the are computed so that the item waves arrive together at time periods 
1, . . . , n — 1 and advance linearly inbetween. 
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Algorithm JRPPD{n,K,P € Z+; /° g M+; /,c e Rf^; 

(1) Start with the wave a,t W = n and the dual solution {v,w,u) = 0. All 
orders are closed, and all demand points are unserved, i.e. {x, y) = 0. 

(2) While there are unserved demand points: 

(3) Decrease W continuously and update W'^ according to (36). At the same 
time, for unserved demand points (<, fc), increase v'l = max{0, dj^(/iit — 
W''), and update Wg-nt and Usi^t so as to maintain (35). The wave stops 
when a joint or individual order becomes tight. 

(4) If an individual order (s,fc,7r) became tight, fix the variables w'l^^ as 
described in (35). If the joint order (s, tt) is also tight, serve all demand 
points contributing to (s, fc,7r) from (s,7r). 

(5) If a joint order (s, tt) became tight, open the joint order and all tight 
individual orders (s.n.k). For each unserved demand point {t,k) that 
contributes to joint order (s,7r), serve {t,k) from (s,7r). 

(6) For each open joint order s from 1 to n: 

(7) If there is a demand point (f, k) that contributes to s and to another 
open joint order s' with s' < s, close s. 

(8) For each item k: 

(9) While not all demand points have been processed in step (11): 

(10) Select the latest such demand point {t,k). Let ireeze{t,k) be the 
location of W'^ when was stopped, and let s be the earliest open 
joint order in [freeze(t, fc), i]. 

(11) Open individual order {s,k). Serve all demand points {t',k) with 
s <t' <t from (s, k). 

(12) Return (a;, y) and (w,w). 

A direct implementation of this algorithm will have an exponential running time. It 
is possible to implement this algorithm to have a polynomial running time, however we 
will not do so here. Instead, we only prove that it provides a 4-approximation guarantee. 

The proof closely resembles that from [LRS06] , and therefore for this proof we assume the 
reader is fully familiar with the joint replenishment results from [LRS06]. 

Theorem 5. Algorithm JRPPD provides a ^-aVPfoximaUon guarantee for the generalized 
JRP. 

Proof. First, similarly to the proof of Levi et al. and Theorem 3, after loop (6), each 
demand point contributes to at most one open joint order. Since we do not open any 
other joint orders after this step, the joint order cost is fully paid by the dual solution, i.e. 
Z^"=i '}2i-K(^[p]K f'jJsTT — SfcLi SiLi ^t^- °f ^ times the cost of the dual solution, we 

allocate one toward the cost of the joint orders. Therefore, we need not consider the cost 
of the joint orders further in this proof. 

Second, also similarly to the proof of Levi et al. and Theorem 3, after loop (6), for each 
demand point {t, k) there is at least one open joint order in [freeze(t, k),t\. Therefore, after 
loop (8), the algorithm produces a feasible primal solution. 

Since we have already covered the cost of joint orders, we now consider each item k 
separately. We bound the holding cost and the cost of individual orders in terms of the dual 
value, similarly to Levi et al. Due to the different cost structure of the JRP and generalized 
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JRP, we are only able to bound the holding and individual order cost by 3 times the cost 
of the dual solution, i.e. ^^=1 E.e[P]- f^Vs. + E:=i Ilts E.e[P]-(4 + h'^^x'^^ < 

Therefore, we obtain a 4 approximation algorithm. □ 
4.2 Applying the Technique 

Finally, we obtain the strongly polynomial algorithm for the concave cost JRP. First, we 
reduce the concave cost JRP to an infinitely-sized generalized JRP: 

E + E fLVs. + E (4 + (37a) 



mm 



s6[n] se[n],ke[K] l<s<t<n 



s.t. E ^'^t = 1' l<t<n,ke [K], (37b) 
se[t] 

0<xJ^t<yL, l<s<t<n,ke[K],TreR'l, (37c) 

< < y,^,, l<s<t<n,ke[KUe M:^, (37d) 

G {0, 1}, y,^, G {0, 1}, 1 < s < n,k e [K],7r e . (37e) 

As before, let ConcaveJRPPD be the algorithm obtained by executing Algorithm JRPPD 
implicitly on infinitely-sized problem (37). 

Theorem 6. Algorithm ConcaveJRPPD is a strongly polynomial 4- approximation algo- 
rithm for the concave cost JRP. 

Proof. Although in this setting all ordering costs are the same over time, we will need to 
refer to ordering costs and groups of tangents at specific times. With this in mind, we will 
refer to the ordering cost of item k at time t by {(j)^ ,t) and to the joint ordering cost at 
time t by (/,t). 

Note that we do not need to keep track of variables for each tt G explicitly. Denote 
the tangent to the individual ordering cost {(f)^ , s) that becomes tight first by pj^^. Then all 
the other tangents to this individual ordering cost at this time are no longer relevant: 

1. Concerning individual ordering costs. For any item I ^ k, the behavior of demand 
points {t, I) or tangents to costs {(j)\t) does not depend on item k, except through the 
joint ordering cost. 

2. Concerning the joint ordering cost. For any wave position, the contribution to the 
joint ordering cost EfeLi Ylt=s '"'sTvt highest for tt with = p^^. 

Therefore, it suffices to keep track, for each item k and time s, of the wave position when 
the first tangent to {(f)^,s) becomes tight. When this occurs, we can stop considering all 
other tangents to {(f)^,s). When computing the wave position when the joint ordering cost 
becomes tight, we need to consider only the tangents that became tight for individual 
ordering costs (^*^,s). Through this transformation, the wave position when the joint 
ordering cost becomes tight can be computed by Lemma 6. 

We now define the following events and wave positions when they occurred: 
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Wave Pos. Event 



Wi{t) The wave reaches time period t, i.e. W = t. 

W2{t, k) A tangent p of order point (<, k) becomes tight. 

W-i{t) The joint order at time t becomes tight. 

The computation now proceeds similarly to the lot-sizing case. We compute the largest of 
the wave positions W\{t), W2(t, fc), and Wz{t) (which corresponds to the smallest time in 
the facility location problem). After the computation we update the other Vl^- values, and 
iterate. □ 
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